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Abstract. Results will be presented from a study of the QCD transition region using 2+1 flavors of fermions and a dislocation 
suppressing gauge action on a lattice with temporal extent of 8 and spatial extent 16 (1.9 — 2.7 fm). A series of temperatures 
from 140 through 200 MeV, separated by 10 MeV have been studied. All the simulations lie on a line of constant physics 
with 200 MeV pions, realized using domain wall fermion, a chirally symmetric fermion formulation. The chiral condensates, 
susceptibility, anomalous symmetry breaking and a detailed study of the Dirac spectrum will be described and compared with 
earlier staggered results. 
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The chiral phase transition of strongly-interacting matter has been an intriguing topic for decades. Many significant 
efforts, both theoretical and experimental, have been dedicated to understanding this phenomena. Lattice Quantum 
Chromodynamics, since the day of its birth, has offered a powerful technique for an ab initio, non-perturbative study 
of the transition region. Here we present a preliminary study of the transition region using a chiral formulation of 
lattice fermions in order to more accurately realize the chiral symmetry responsible for the transition. 

GAUGE AND FERMION ACTIONS 

While the gauge field can be implemented on the lattice in a straightforward manner, the usual lattice fermion 
schemes suffer from problems such as breaking of chiral symmetry or creating doublers etc., which interfere with the 
exploration of the QCD chiral phase transition using lattice methods. The domain wall fermion (DWF) formulation 
[1, 2], a variant of Wilson fermions, avoids these problems at the cost of introducing an extra, s dimension. The left 
and right-handed chiral states are bound to the four dimensional boundaries (walls) at the two ends of the 5-D volume. 
Unphysical, massive degrees of freedom from the fifth dimension are to a large extent canceled by the introduction of 
a Pauli-Villars term in the simulation. The low modes produce a discrete, four-dimensional version of QCD. 

The DWF formulation has chiral breaking symmetry under good control. The residual chiral symmetry breaking 
is characterized by a residual, additive quark mass: OTres(ii) = ciexp(— AcAs)/As +C2/Ls- If the extension in fifth 
dimension, L, is taken to infinity, one will obtain the overlap fermion formulation which exactly preserves chiral 
symmetry. The infinite Ls limit commutes with the continuum extrapolation, allowing us to choose an appropriate 
Ls to circumvent the heavy computational burden of overlap fermions while keeping the chiral symmetry breaking 
minimal. Another benefit of the DWF formalism is that U{1)a symmetry is only broken by axial anomaly in contrast 
to staggered fermions, where the anomalous symmetry is also broken by a^ artifacts. 

In the simulation of QCD thermodynamics on a lattice with constant physical quark masses, one naturally needs a 
small bare quark mass. Since the residual mass is an additive correction to the bare quark mass in DWF formalism, 
viz. mq = mi-es +»^input, we must also minimize mres- From expression for m^es, the residual mass can be decomposed 
into two parts. The first part is produced by standard five dimensional states and is exponentially suppressed; the other, 
comes from gauge field dislocations associated with changing topology and is only suppressed by l/Ly. Therefore, 
naively increasing will not only be computationally expensive but also inefficient. Fortunately, the topology change 
associated with the second term also leads to zero modes of an unphysical 4-dimensional Dirac operator, M5) 
so we could multiply a factor (e.g. determinant of that Dirac operator) to suppress those configurations. Since such a 



factor will freeze topological tunneling, we add a twisted mass term to the determinant and introduce the ratio; 
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This is the dislocation suppression determinant ratio (DSDR) [3, 4, 5]. This alters predominately the ultraviolet part of 
the theory, so the physical quantities in which we are interested are minimally affected. 

TABLE 1. Summary of the ensembles. 
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Table 1 summarizes the ensembles in our simulation. We use 2+1 flavors of DWF with Iwasaki and DSDR gauge 
action. The spatial extent of the lattice is 16 (about 1.9 - 2.7 fm). The temporal extent is fixed to be 8 and the j3 values 
are chosen so that we have a range of temperatures from 140 to 200 MeV. The bare quark masses are also adjusted to 
make our ensembles lie on a line of constant physics with a pion mass of about 200 MeV. 



CHIRAL SYMMETRY 

The first quantity of interest is the chiral condensate shown as a function of Monte Carlo time in fig. 1 . As compared 
to those below 160MeV, the evolutions of the chiral condensate for temperatures above 160MeV display a distinct 
character, fluctuating above a base line. The graph of the disconnected susceptibilities gives a more direct indication 
of a transition temperature around 160 MeV, which agrees quite well with the staggered results. However, since our 
spatial volume is quite limited (an aspect ratio 2), a quantitative comparison is likely premature. 




Trajectory T (McV) 

FIGURE 1. Evolution of {ijnif)i, each shifted up by 12 units (left) and preliminary results of the disconnected susceptibility, 
normalized in a consistent scheme differing from that used in Table 1, at various temperatures (right). 

The chiral condensate is closely related to the small eigenvalues of the Dirac operator, which thus provides an 
alternative perspective into the symmetries associated with the QCD phase transition. For instance, the well-known 
Banks-Casher relation [6], lim„j,^olimv^oo(v//V^/) = — 7rlim;L^olinim/^olimv^ooP(A), relates the chiral condensate 
and the A = value of the density p(A) of Dirac eigenvalues. 



The Dirac spectrum also provides unique insights into the anomalous U{1)a symmetry. Specifically, the U{1)a- 
breaking difference of pseudoscalar and scalar susceptibilities is given by: 



(2) 



While Ajj;_s is expected to receive a small, dilute instanton gas contribution from p{X)^ 5 (A ) at large temperatures, 
there may also be larger, lower-temperature non-semi-classical contribution from p(A) cx A. We used Kalkreuter- 
Simma's method [7] to calculate the lowest 100 eigenvalues of the hermitian version of DWF Dirac operator. The 
eigen-spectrum is then renormalized using an method similar to Giusti and Luscher [8]. 
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FIGURE 2. Dirac spectrum of T = 150 - 200 MeV ensembles. 

As can be seen in fig. 2, the renormalization aligns the light quark mass (dashed vertical line) with the near-zero 
mode peak. It also compensates for the discrepancy in the Banks-Casher relation, which however still does not agree 
at about 50% level below the transition. We attribute the discrepancy to the finite volume or finite mass effects. Above 
170MeV, the chiral condensate and eigendensity both start to vanish as expected for temperatures above the transition. 
But it remains unclear whether the slope of the eigendensity vanishes at ISOMeV. At even higher temperature, a 
possible gap starts to emerge, indicating an effective restoration of f/(l)A symmetry. 

We conclude that DWF provides a well-suited tool to study the QCD phase transition. From the chiral susceptibility, 
we find indication of a transition temperature around 160MeV. Well above the transition, we observe possible signals 
of t/(l)A restoration. Larger volume results and chiral extrapolation will be needed for a more definite conclusion. 
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